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Abstract 



The first calculation of the gluon propagator using an 0{a 2 ) improved ac- 
tion with the corresponding 0{a?) improved Landau gauge fixing condition 
is presented. The gluon propagator obtained from the improved action and 
improved Landau gauge condition is compared with earlier unimproved re- 
sults on similar physical lattice volumes of 3.2 3 x 6.4 fm 4 . We find agreement 
between the improved propagator calculated on a coarse lattice with lattice 
spacing a = 0.35 fm and the unimproved propagator calculated on a fine 
lattice with spacing a = 0.10 fm. This motivates us to calculate the gluon 
propagator on a coarse large-volume lattice 5.6 3 x 11.2 fm 4 . The infrared 
behavior of previous studies is confirmed in this work. The gluon propagator 
is enhanced at intermediate momenta and suppressed at infrared momenta. 
Therefore the observed infrared suppression of the Landau gauge gluon prop- 
agator is not a finite volume effect. 
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I. INTRODUCTION 



There has been considerable interest in the infrared behavior of the gluon propagator, 
as a probe into the mechanism of confinement and as input for other calculations. Lattice 
gauge theory is an excellent means to study such nonperturbative behavior. See, for example, 
Ref. [IJ for a recent survey. 

The infrared part of any lattice calculation may be affected by the finite volume of the 
lattice. Larger volumes mean either more lattice points (with increased computational cost) 
or coarser lattices (with corresponding discretization errors). The desire for larger physical 
volumes thus provides strong motivation for using improved actions. Improved actions have 
been shown to reduce discretization effects 0, although some concerns have been expressed 
that coarse lattices may miss important instanton physics. 

In this study, no change is seen in the infrared gluon propagator, even with a lattice 
spacing as coarse as 0.35 fm. We find the gluon propagator to be less singular than in 
the infrared. Our results suggest that the gluon propagator is infrared finite, although more 
data is needed in the far infrared to be conclusive. This behavior is similar to that observed 
in three-dimensional SU(2) studies M. 

II. 0(a 2 ) IMPROVEMENT 

The 0(a 2 ) tadpole- improved action is defined as 

s g = ^rJ2 neTr ( 1 ~ U A*)) - T^E^ eT Kl - U iect (x)), (2.1) 

6 pi llU rect 

where the operators U p \(x) and U TCCt (x) are defined 

U pl (x) = U^x)U v (x + jj)Ul(x + u)Ul(x), (2.2) 

and 

^rectO) = U^(x)U u (x + (l)U v (x + i> + foU^X + 2fi)Ul(x + v)Ul(x) 

+ U^(x)U^x + fi)U v {x + 2p)Ul(x + fi + v)Ul(x + v)Ul{x). (2.3) 

The link product U iect (x) denotes the rectangular 1x2 and 2x1 plaquettes. For the tadpole 
(mean-field) improvement parameter we use the plaquette measure 

u = {hleTr < U pl >) * . (2.4) 

Eq. ( |2.1D reproduces the continuum action as a — > 0, provided that (3 takes the standard 
value of 6/g 2 . Note that our (3 = 6/g 2 differs from that used in PH[7J. Multiplication of 
our (3 in Eq. ( |2.1|) by a factor of 5/3 reproduces their definition. 0(g 2 a 2 ) corrections to this 
action are estimated to be of the order of two to three percent M. 

Gauge configurations are generated using the Cabbibo-Marinari || pseudoheat-bath al- 
gorithm with three diagonal SU{2) subgroups. The mean link, u , is averaged every 10 
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sweeps and updated during thermalization. Representative gauge fields are selected after 
5000 thermalization sweeps. 

Gauge fixing on the lattice is achieved by maximizing a functional, the extremum of 
which implies the gauge fixing condition. The usual Landau gauge fixing functional || is 

^[{U}} = £ ^Tr {U°(x) + U°(xy} , (2.5) 



where 



and 



U%(x) = G{x)U IJ ,{x)G{x + fi)\ (2.6) 



G(x) = exp j-iX> a (:r)T a j. (2.7) 

A maximum of Eq. ( |2.5| ) implies that J^/id^A^ = up to 0(a 2 ). To ensure that gauge 
dependent quantities are also 0(a 2 ) improved, we implement the analogous 0(a 2 ) improved 
gauge fixing functional 



as described in |I(J. We employ a Conjugate Gradient, Fourier Accelerated gauge fixing 
algorithm [|ll|] optimally designed for parallel machines. 

III. THE LANDAU GAUGE GLUON PROPAGATOR 

The gauge links U^x) are expressed in terms of the continuum gluon fields as 

U^x) =Ve i90 ^ A ^ z)dz , (3.1) 

where V denotes path ordering. From this, the dimensionless lattice gluon field A^(x) may 
be obtained via 

M* + A/ 2 ) = ^ fax) - Ul(x)) - ^-Tr (u,(x) - U^x)) , (3.2) 



accurate to 0(a 2 ). 0(a 2 ) improved gluon field operators have also been investigated. While 
the infrared behavior is unaffected by the improvement, the ultraviolet behavior suffers due 
to the extended nature of an improved operator. These results will be discussed in detail 
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We calculate the gluon propagator in coordinate space 

Df u (x,y) = (A;(x)Al(y)), (3.3) 
using (pO). To improve statistics, we use translational invariance and calculate 



K(y)-U^ A > )A ^ x+y)) - (3 - 4) 

V X 

In this report we focus on the scalar part of the propagator, 

D (y) = tv^t^ mziT.DZiv)' ( 3 - 5 ) 

where Nd = 4 and N c = 3 are the number of dimensions and colors. This is then Fourier 
transformed into momentum space 



£(?) = £e*»£>(y), (3.6) 

y 

where the available momentum values, q, are given by 

= "><- L t L t\- ^ 

and L M is the length of the box in the \x direction. In the continuum, the scalar function 
D(q 2 ) is related to the Landau gauge gluon propagator via 

DfM = - qj fW b Dtf) . (3.8) 

The bare, dimensionless lattice gluon propagator, D(qa) is related to the renormalized 
continuum propagator, Dn(q; ji) via 

a 2 D(qa) = Z 3 (fj,, a)D R (q; /i), (3.9) 

where the renormalization constant, Z^(fi,a) is determined by imposing a renormalization 
condition at some chosen renormalization scale, //, e.g., 

D R (q)\ q2= ^ = — 2 . (3.10) 

This means that there is an undetermined multiplicative renormalization factor, Z^(a), on 
each of our propagators. Since our purpose is to compare our two improved, coarse lattices 
with the unimproved, finer one, it is sufficient to consider only the their relative renormal- 
izations. We have slightly rescaled the improved propagators so as to provide a reasonable 
match with old one. The relevant quantity is 

Z 3 (0.10)/Z 3 (0.35) = 1.02. (3.11) 
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IV. ANALYSIS OF LATTICE ARTIFACTS 



To emphasize the nonperturbative behavior of the gluon propagator, we divide the prop- 
agator by the tree level result of lattice perturbation theory. Hence Figs. [1], and |3] are 
plotted with q 2 D(q 2 ) on the y-axis, which is expected to approach a constant up to loga- 
rithmic corrections as q 2 — > oo. Here 

? M =-sm ~y . (4.1) 

To reduce ultraviolet noise resulting from the lattice discretization, the available mo- 
menta are cut half way into the Brillouin zone, that is 

gmax = ■ (4.2) 

2a 

All figures have a cylinder cut imposed upon them, i.e. all momenta must lie within a cylinder 
of radius two spatial momentum units centered about the lattice diagonal. The propagators 



are plotted in physical units, which we obtain from the string tension [13| with yfo = 440 
MeV. 

Fig. |I] is reproduced from Ref. [1~4[, where the standard, Wilson action is used. The 
propagator is calculated on a 32 3 x 64 lattice at j3 — 6.0, which corresponds to a lattice 
spacing of 0.1 fm. This propagator produces the correct asymptotic behavior and a detailed 
analysis shows the anisotropic finite volume errors are small. However, it was impossible to 
rule out isotropic finite volume artifacts. 

We use the improved action described above to calculate the gluon propagator on a small 
(10 3 x 20), coarse (a = 0.35 fm) lattice, which is shown in Fig. Q Despite the coarse lattice 
spacing we see that it reproduces the infrared behavior of Fig. [I]. 

Finally, we calculate the propagator on a 16 3 x 32 lattice, at the same (3 providing 
a = 0.35 fm. This corresponds to a very large physical volume of 5.6 3 x 11.2 fm 4 . Fig. |3| 
illustrates the results. These results largely agree with the previous two calculations of the 
propagators. The behavior of the gluon propagator is not changed by changing the volume. 

In Fig. ^ we have superimposed the gluon propagators for all three lattices. Here we 
plot D(q 2 ) on the y-axis to allow an alternative examination of the most infrared momenta. 
The points at > 350 MeV are very robust with respect to volume. Only the very lowest 
momenta points show signs of finite volume effects. With more volumes it should be possible 
to extrapolate to the infinite volume limit. We see that the propagator, at the very lowest 
momentum points, decreases as the volume increases. This strongly suggests an infrared 
finite propagator. 



V. CONCLUSIONS 

The gluon propagator has been calculated on a coarse lattice with an 0(a 2 ) improved 
action, in 0(a 2 ) improved Landau gauge. The infrared behavior of this propagator is con- 
sistent with that of a previous study on a finer lattice with an unimproved action, but 
comparable volume. The propagator was then calculated on another improved lattice with 
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FIG. 1. Gluon propagator from 75 standard, Wilson configurations, on a 32 3 x 64 lattice with 
spacing a = 0.10 fm. 
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FIG. 2. Gluon propagator from 75 tree-level improved configurations on a 10 3 x 20 lattice with 
spacing a = 0.35 fm, and a physical volume of 3.5 3 x 7 fm 4 . 
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FIG. 3. Gluon propagator from 75 tree-level improved configurations, on a 16 3 x 32 lattice with 
spacing a = 0.35 fm, and a physical volume of 5.6 3 x 11.2 fm 4 . 
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FIG. 4. Comparison of the gluon propagator on the three different lattices. The volumes are 
3.2 3 x 6.4 fm 4 , 3.5 3 x 7.0 fm 4 , and 5.6 3 x 11.2 fm 4 . 
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the same spacing, but larger volume. The increase in volume left the propagator largely 
unchanged. In particular, it has been shown that the turnover observed in 14] is not a finite 
volume effect. 

With more lattices it may be possible to extrapolate to infinite volume, but from this 
study we can only make tentative conclusions. We have ruled out the q~ 4 behavior popular 
in Dyson- Schwinger studies, and any infrared singularity appears to be unlikely. An infrared 
finite propagator is most plausible. The gluon propagator would need to drop rapidly for 
momenta below ~ 350 MeV in order to vanish as suggested by Zwanziger |15j and others. 
Even larger volume lattices will be needed to study this possibility. The possible effects of 
lattice Gribov copies remains a very interesting question and we plan to study this in the 
near future. 
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